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Abstract; In this letter we demonstrate that our dipole model is successful in describing the inclusive 
production within the same framework as diffractive physics. We believe that this achievement stems from 
the fact that our approach incorporates the positive features of the Reggeon approach and CGC/saturation 
effective theory, for high energy QCD. 
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1 Introduction 

The LHC data on inclusive production [1-3] call for a theoretical understanding of these processes within the 
framework of QCD. At first sight it would appear, that this process is a typical soft process, which occurs at 
long distances, where one should use the methods of non-perturbative QCD. Since such methods are only in 
an embryonic stage, soft processes at high energy remain in the arena of high energy phenomenology, based 
on the concept of a soft Pomeron. Adopting this approach, inclusive production can be calculated using 
the technique of Mueller diagrams [4], It has been demonstrated that soft Pomeron based models, provide 
a reasonable decription of the data[5, 6]. The advantage of our approach, is the feasibility of describing 
inclusive production on the same footing as diffractive production, and elastic scattering. 

On the other hand, in the CGC/saturation approach for iclusive production[7-13], one has a different 
senario. In this approach the inclusive production occurs in two stages. The hrst stage is the production of 
a mini-jet with the typical transverse momentum Qg, where Qs is the saturation scale, which is much larger 
than the soft scale. This process is under full theoretical control. The second stage is the decay of the mini¬ 
jet into hadrons, which has to be treated phenomenologically, using data from the hard processes. Such 
approach leads to a good description of the experimental data on inclusive production, both for hadron- 
hadron, hadron-nucleus and nucleus-nucleus collisions, and observation of regularities in the data, such as 
geometric scaling [14-18]. The shortcoming of this approach is the fact that it is detached from diffractive 
physics. 

It should be mentioned, that the recently published measurements of the pseudorapidity distributions 
of charged particles in proton-proton collisions at an energy of 8 TeV, provide an additional challenge for 
model builders, which has not yet been successfully answered [2]. 

In this letter, we continue (see Refs. [19, 20] to construct a model for high energy soft interactions, 
which incorporates the advantages of both approaches. This model is based on the Colour Glass Con- 
densate(CGC)/saturation effective theory (see Ref.[21] for the review), and on the perturbative BFKL 
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Pomeron[22]. We assume that the unknown mechanism for the confinement of quarks and gluons in QCD, 
is not important, and its influence can be reduced to the determination of several parameters related to the 
CGC/saturation approach, which depend on long distance physics. 

The main attributes of the model have been discussed in Refs. [19, 20], in this paper we will only include 
information that we require for the discussion of inclusive production. 

2 Main formulae 

First, we discuss the initial stage of hadron production in the framework of CGC/saturation approach. For 
mini-jet production, we use the kx factorization formula, that has been proven in Ref. [23] ( see aslso Refs. 
[24-29] where this proof has been verified). 

= ^ / “"•‘T 4>'3 (x#T - St) (2.1) 

where (f)^ denotes the probability to find a gluon that carries the fraction Xi of energy with k_i_ transverse 
momentum, and as = asNc/ir, with the number of colours equal to W- + y = ln(l/xi) and 
— y = ln(l/x 2 ). (j)Q is the solution of the Balitsky-Kovchegov(BK) [10, 12] non-linear evolution equation, 
and can be viewed as the sum of ‘fan’ diagrams of the BFKL Pomeron interactions, shown in Fig. 1. 



Figure 1. The graphic representation of Eq. (2.1) (see Fig. l-a).For the sake of simplicity all other indices in 
cj) {xi,pt — kx) and (j>{x 2 ,kT) are omitted. The wavy lines denote the BFKL Pomerons, while the helical lines 
illustrate the gluons. In Fig. 1-b the Mueller diagram for the inclusive production is shown. 

In our model we use the simple formula which is a good approximation to the numerical solution of 
the BK equation, see Ref. [30]: viz. 

[Gjp {z)) = a (1 - exp {-Gjp {z))) + (1 - a) . ' , (2.2) 

1 -I- Gjp [z) 
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with a = 0.65 and G]p{z) = 4 >q where we have used two inputs: r = R and Qg = 

S {mb) exp (Ay). For the values of 1 — jcr and A, we have estimates in the leading order of 
perturbative QCD: 1 — 'jcr = 0.63 and A = 4.880:5. The value of A is a fitting parameter, which effectively 
includes the higher order QCD corrections. In this paper we use the value A = 0.38 which we found in 
Ref. [20]. The parameter m and the function S (mb) originate from non-perturbative QCD contributions, 
and are taken in the following form: 

2 p 

S {mb) = ^6“™^ where / cPbS{b) = 1 (2.3) 


(pQ can be calculated using the initial conditions, of the BFKL equation. 

However, we do not know these conditions, and so we consider (/>o as an additional phenomenological 
parameter. The values of these parameters are taken from Ref.[20]: m = 5.25GeR and (pQ = 0.0019. All 
these parameters describe the CGC /saturation structure of the BFKL Pomerons and their interactions. 
We introduce phenomenological parameters to describe the structure of the hadron. We choose the two 
channel model for such a structure and describe the vertex of the BFKL Pomeron interaction with a hadron 
state i in the following form: 


gi{mi,b) = giSip{mi,b) where Sip{mi,b) 


dvr 


b Ki {mi b) 


(2.4) 


The parameters that we use in this paper, have been extracted from fitting the elastic and diffractive data 
in Ref. [20], and their values are: 


= 110.2mi = 0.92GeK; ^ 11.2 m2 = 1.9GeK; (2.5) 


Finally, Eq. (2.1) can be re-written as a Mueller diagram of Fig. 1-b, and the inclusive cross section is 
given by 


^ = arrln(H'/ir„)|a-‘/„«(iF+!,) /nWQr-s) 

-h Qy ~ ) 

+ ( 3 ^ Qy + y) /n(2) Qy - y^ | (2.6) 

where a and /3 describe the structure of the diffractive scattering in the two channel model, where the 
observed physical hadronic and diffractive states are written in the form 

iph = Q:'I'i + ;d^2; 'pD = —/3'Ll + a ^^2; where +/3^ = 1; (2-7) 

Inb) is given by 

(y) = I d% 5 (m^, b) Gjp (y)) (2.8) 
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where G]p{y) = (/)oexp(A(l —Jcr)y) and is defined in Eq. (2.2). Regarding the factor in front 

of Eq. (2.6) i.e. \n.{W/W q), where W = y/s is the energy of collision in c.m.f., and Wq is the value of 
energy from which we can start our approach. One can see that Eq. (2.1) is divergent in the region of 
small pt < Qs- Indeed, in this region (/>’s in Eq. (2.1) do not depend on px, since kx ^ Qs > PT, and 
the integration over px leads to In where rusoft is the non-perturbative scale that includes the 

confinement of quarks and gluons {irisoft ~ ^qcd)- 

To convert the rapidity distribution (which we calculate theoretically), to pseudorapidity ( 77 ) one, we 
need to know the mass of min-jet (m). The simple formulae are well known (see Ref. [14] for example): 


y{ri,PT) = 2^^' 


—- + smh r] + sm 77 
Pt 


+ sinh^ 77 — sinh 77 


(2.9) 


with the Jacobian 


h{v,px) = 


cosh 77 




-|- sinh^ 77 


( 2 . 10 ) 


The mass of mini jet is given by = 2msoftPT (see Ref. [14]). Since the typical transverse momentum 
is equal to the saturation scale, we have 


_ 2msoft _ 2 (^\ 

pI Qs{W) "'^WoJ 


( 2 . 11 ) 


where Tq and ojpjp are phenomenological parameters that are determined from the experimental data. 


Finally, 

^ = h{r],Qs)^{y{v,Qs)) 


( 2 . 12 ) 


3 Comparison with the experimental data 

In Fig. 2 we plot our predictions compared to the experimental data. As we have mentioned all other 
parameters have been extracted from the diffractive and elastic data in Ref. [20]. The only free parameters 
are ajpjp and rg. The curves in Fig. 2 are calculated for apjp = 0.21 and rg = 8. 

From Eq. (2.6), we note that the inclusive cross section is sensitive to the contribution of the black 
component. As we discussed in Ref. [20], qualitatively, we have in our two channel model two different 
components: one which is transparent, even at ultra high energy (e.g. at IF = 57TeV) while the second 
component, starts being black at rather low energy (say at IF = 0.9 TeF). Hence, our good description of 
the experimental data, checks that the value of this component is consistent with the inclusive data. 

One can see that our model describes the value of the inclusive densities p = {l/ax!SD)d<^/dp and 
their dependence on energy and rapidity, rather well. It should be stressed that the values for gmsd were 
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Figure 2. The single inclusive density ( {I/ctnsd) da/drj) versus energy. The data were taken from Refs. [1-3] and 
from Ref. [31]. description of the multi-particle production in hard processes in The description of the CMS data is 
plotted in Fig. 2-a, while Fig. 2-b presents the comparison with all inclusive spectra with W > 0.546 TeR. 

calculated in our model. Regarding the new data at W= 8 TeV (Fig. 2 and 3), the comparison shows that 
the result of our approach, is slightly below the experimental central values, while the numerous Monte 
Carlo simulations overshoot the data in the central region (see Fig. 3-a). 

Note, we have only dealt with data in the central region, since we do not take into account parton 
correlations due to energy conservation. These are important in the fragmentation region, but difficult to 
incorporate in our present framework. 

In Fig. 3-b it is shown the energy dependence of dNch/dr] at r/ = 0. In the CGC/saturation approach 
dNch /dri\^=Q (X W^, where A corresponds to the energy dependence of the saturation scale. In our model 
the energy dependence is more complicated and can be approximate as Note, that the power 0.29 

is much less than the value of A (A = 0.38). Recall, that in our model A characterizes the energy behavior 
of the saturation scale. 

Concluding this discussion we see that we reached a good description of the data but our estimates are 
below of the data at small values of rapidity r]. We believe, that this is a reflection of our simplified relation 
between y and y. 
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Inclusive pp CMS-TOTEM, is = 8 TeV, L = 45iib ’ 
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2 — >1 in 5.3<ti<6.5 or -6.5<ri<-5.3 



Fig. 3-a 



Fig. 3-b 


Figure 3. The comparison of the inclusive production at VF = 8 TeV with the Monte Carlo models is shown in 
Fig. 3-a. The picture is taken from Ref.[2]. In Fig. 3-b it shown dNch/drj at 77 = 0 versus energy W. Our estimates 
are shown by the solid line. The dotted line corresponds to fit: 0.725 (VF/Wo)° with Wq = 1 GeV. The data are 
taken from Refs.[l, 2, 32, 33]. 


4 Conclusions 


In this letter we demonstrate that our model for the soft (long distance) interaction which is based on 
CGC/saturation approach, is able to describe inclusive production. In other words, we give the example 
that our model can describe both the diffractive (elastic) physics at high energy and the typical produc¬ 
tion process which in the majority of the approaches are treated in different ways. We believe that such 
a successful approach roots in the fact, that our procedure incorporates the advantages of high energy 
phenomenology based on the soft Pomeron interactions , and of the CGC/saturation effective theory that 
includes the description of multi-particle production in perturbative QCD. 

On the other hand, this letter is the natural next step in our search for a model based on QCD, that will 
be able to describe the typical properties of high energy interaction and to include the diffractive production 
and multi-particle generation process on the same footing. It is also a next step in our attempts to build 
this description without addressing the Monte Carlo simulation methods. 


Acknowledgements 

We thank our colleagues at Tel Aviv university and UTFSM for encouraging discussions. Our special thanks 
go to Carlos Contreras, Alex Kovner and Misha Lublinsky for elucidating discussions on the subject of this 


6 















paper. This research was supported by the BSF grant 2012124 and the Fondecyt (Chile) grant 1140842. 


References 

[1] K. Aamodt et al. [ALICE Collaboration], Eur. Phys. J. C 68 (2010) 89 [arXiv: 1004.3034 [hep-ex]]; ALICE 
Collaboration, Eur. Phys. J. C 65 (2010) 111 ]arXiv:0911.5430 ]hep-ex]]. 

]2] S. Chatrchyan et al. [CMS and TOTEM Collaborations], Eur. Phys. J. C 74 (2014) 10, 3053 ]arXiv:1405.0722 
]hep-ex]]; V. Khachatryan et al. [CMS Collaboration], Phys. Rev. Lett. 105, 022002 (2010; ]arXiv:1005.3299 
[hep-ex]]. V. Khachatryan et al. [CMS Collaboration], JHEP 1002 (2010) 041 ]arXiv:1002.0621 ]hep-ex]]. 

]3] ATLAS Collaboration, arXiv: 1003.3124 [hep-ex]. 

]4] A. H. Mueller, Phys. Rev. D2 (1970) 2963. 

]5] E. Gotsman, E. Levin and U. Maor, Phys. Rev. D 81, 051501 (2010) ]arXiv:1001.5157 ]hep-ph]]. 

]6] E. Gotsman, E. Levin and U. Maor, Phys. Rev. D 84, 051502 (2011) [arXiv: 1103.4509 ]hep-ph]]. 

[7] L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. Rep. 100 (1983) 1. 

[8] A. H. Mueller and J. Qiu, Nucl. Phys. B268 (1986) 427. 

[9] L. McLerran and R. Venugopalan, Phys. Rev. D49 (1994) 2233, 3352; D50 (1994) 2225; D53 (1996) 458; 
D59 (1999) 09400. 

[10] 1. Balitsky, ]arXiv:hep-ph/9509348]; Phys. Rev. D60, 014020 (1999) ]arXiv:hep-ph/9812311] 

[11] A. H. Mueller, Nucl. Phys. B415 (1994) 373; B437 (1995) 107. 

[12] Y. V. Kovchegov, Phys. Rev. D60, 034008 (1999), ]arXiv:hep-ph/9901281]. 

[13] J. Jalilian-Marian, A. Kovner, A. Leonidov and H. Weigert, Phys. Rev. D59, 014014 (1999), 
[arXiv:hep-ph/9706377]; Nucl. P/iys. B504, 415 (1997), ]arXiv:hep-ph/9701284]; J. Jalilian-Marian, 

A. Kovner and H. Weigert, Phys. Rev. D59, 014015 (1999), ]arXiv:hep-ph/9709432]; A. Kovner, 

J. G. Milhano and H. Weigert, Phys. Rev. D62, 114005 (2000), ]arXiv:hep-ph/0004014]; E. lancu, 

A. Leonidov and L. D. McLerran, Phys. Lett. B510, 133 (2001); ]arXiv:hep-ph/0102009]; Nucl. Phys. A692, 
583 (2001), ]arXiv:hep-ph/0011241]; E. Ferreiro, E. lancu, A. Leonidov and L. McLerran, Nucl. Phys. A703, 
489 (2002), ]arXiv:hep-ph/0109115]; H. Weigert, Nucl. Phys. A703, 823 (2002), ]arXiv:hep-ph/0004044]. 

[14] D. Kharzeev, E. Levin and M. Nardi, Nucl. Phys. A 730 (2004) 448 [Erratum-ibid. A 743 (2004) 329] 
[arXiv:hep-ph/0212316]; Phys. Rev. G 71 (2005) 054903 ]arXiv:hep-ph/0111315]; D. Kharzeev and E. Levin, 
Phys. Lett. B 523 (2001) 79 ]arXiv:nucl-th/0108006]; D. Kharzeev and M. Nardi, Phys. Lett. B 507 (2001) 
121 ]arXiv:nucl-th/0012025]. 

[15] D. Kharzeev, E. Levin and M. Nardi, Nucl. Phys. A 747 (2005) 609 ]arXiv:hep-ph/0408050]. 

[16] E. Levin and A. H. Rezaeian, AIP Conf. Proc. 1350 (2011) 243 ]arXiv:1011.3591 ]hep-ph]]; Phys. Rev. D 82 
(2010) 014022 [arXiv: 1005.0631 ]hep-ph]]; [arXiv: 1102.2385 ]hep-ph]]; Phys. Rev D82 (2010) 054003. 
[arXiv:1007.2430 ]hep-ph]]. 

[17] L. McLerran, M. Praszalowicz, Nucl. Phys. A 916 (2013) 210 ]arXiv:1306.2350 ]hep-ph]]; Acta Phys. Polon. 
B42 (2011) 99 ]arXiv:1011.3403 ]hep-ph]]; Acta Phys. Polon. B41 (2010) 1917. [arXiv: 1006.4293 ]hep-ph]]. 


7 


[18] M. Praszalowicz, Phys. Lett. B 727 (2013) 461 [arXiv: 1308.5911 [hep-ph]]; Acta Phys. Polon. Supp. 6 (2013) 
3, 809 [arXiv:1304.1867 [hep-ph]]; Phys. Lett. B 704 (2011) 566 ]arXiv:1101.6012 ]hep-ph]]l Phys. Rev. Lett. 
106 (2011) 142002 [arXiv: 1101.0585 [hep-ph]]. 

[19] E. Gotsman, E. Levin and U. Maor, Eur. Phys. J. C 75 (2015) 1, 18 [arXiv: 1408.3811 [hep-ph]]. 

[20] E. Gotsman, E. Levin and U. Maor, arXiv:1502.05202 [hep-ph]. 

[21] Yuri V Kovchegov and Eugene Levin, “ Quantum Choromodynamics at High Energies", Cambridge 
Monographs on Particle Physics, Nuclear Physics and Cosmology, Cambridge University Press, 2012 . 

[22] E. A. Kuraev, L. N. Lipatov, and F. S. Fadin, Sou. Phys. JETP 45, 199 (1977); Ya. Ya. Balitsky and L. N. 
Lipatov, Sov. J. Nucl. Phys. 28, 22 (1978). 

[23] Y. V. Kovchegov and K. Tuchin, Phys. Rev. D 65 (2002) 074026 [arXiv:hep-ph/0111362]. 

[24] M. A. Braun, Eur. Phys. J. C 48 (2006) 501 [arXiv:hep-ph/0603060]; Phys. Lett. B 483 (2000), 105. 

[25] J. Jalilian-Marian and Y. V. Kovchegov, Phys. Rev. D 70 (2004) 114017 [Erratum-ibid. D 71 (2005) 079901] 
[arXiv:hep-ph/0405266]. 

[26] C. Marquet, Nucl. Phys. B 705 (2005) 319 [arXiv:hep-ph/0409023]. 

[27] A. Kovner and M. Lublinsky, JHEP 0611 (2006) 083 [arXiv:hep-ph/0609227]. 

[28] E. Levin and A. Prygarin, Phys. Rev. C 78 (2008) 065202 [arXiv:0804.4747 [hep-ph]]. 

[29] A. Kormilitzin, E. Levin and A. Prygarin, Nucl. Phys. A 813 (2008) 1 [arXiv:0807.3413 [hep-ph]]. 

[30] E. Levin, JHEP 1311 (2013) 039 [arXiv: 1308.5052 [hep-ph][. 

[31] C. Amsler et al. (Particle Data Group), Phys. Lett. B667 (2008) 1. 

[32] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. D 79 (2009) 112005 [Erratum-ibid. D 82 (2010) 119903] 
[arXiv:0904.1098 [hep-ex]]; F. Abe et al. [CDF Collaboration], Phys. Rev. D 41 (1990) 2330. 

[33] G. J. Alner et al. [UA5 Collaboration], Z. Phys. C 33 (1986) 1. 


